Dominant Spin-Flip Effects for the Hadronic Produced J/ifj 

Polarization at TEVATRON 

Xing-Gang Wu * and Zhen-Yun Fang 
Department of Physics, Chongqing University, Chongqing 400044; P-R- China 

Abstract 

Dominant spin-flip effects for tlie direct and prompt J /tp polarizations at TEVATRON run II 
with collision energy 1.96 TeV and rapidity cut \y^^'^\ < 0.6, have been systematically studied, 
especially, the spin-flip effect for the transition of (cc)8p5'i] into J/tp has been carefully discussed. 
It is found that the spin-flip effect shall always dilute the J/^|J polarization, and with a suitable 
choice of the parameters ao,i and co,i,2) the J/ip polarization puzzle can be solved to a certain 
degree. At large transverse momentum pt, a for the prompt J/ip is reduced by ~ 50% for /o = 
and by ~ 80% for /o = 1. We also study the indirect J/ip polarization from the 6-decays, which 
however is slightly affected by the same spin-flip effect and then shall provide a better platform to 
determine the color-octet matrix elements. 

PACS numbers: 12.38.Bx, 12.39.Jh, 14.40.Lb 



* email:wuxg@itp. ac.cn 



1 



I. INTRODUCTION 

Within the non-relativistic QCD (NRQCD) framework the hadronic production of 
J/ ip is dominated by the gluon fragmentation in which a gluon fragments into a color-octet 
state (cc)8p5'i]. And if the spin-symmetry hold for charm quarks, as is usually adopted 
in the literature, then the formed J/ip shall always show large transverse polarization at 
sufficiently large pt- But this prediction contradicts with the measured at TEVATRON 
{2, 3|. This is the well-known J/ip polarization puzzle. Recently, the next-to-leading order 
(NLO) analysis of the J/ip polarization have been done by Refs.jj, [sl, with (cc) pair in J/i/j 
Fock expansion being in color-singlet state pS*!]!, color-octets pS'ijg and [^So]s respectively. 
It is found that even with those NLO corrections, the J/ip polarization puzzle still can not 
be solved. This implies that by doing the NLO calculations, one can make the perturbative 
QCD results more convergent on and then more reliable, but one can not change the fact 
that the produced J/tp is largely transversely polarized. 

There are many suggestions to solve such polarization puzzle, e.g. the effects of the 
initial gluon off-shellness may provide a possible explanation for such J/tp spin alignment 
in pp collision Gj, which depends heavily on the unintegrated parton distribution of the 
initial gluon(s). While the spin-flip interaction provides another effective way to solve the 
J/ip polarization puzzle as suggested by Ref. ^. By taking the spin-flip interactions into 
account, it has been argued in Ref. 7|] that the direct J/ip and ip' polarizations can be diluted 
to a certain degree. And an analysis for the prompt ip' production with proper ranges for the 
newly introduced parameters ao,i and Ci^2 has been presented there. In the present paper, we 
shall make a systematical discussion on how the spin-flip interaction affects the polarization 
of the direct, the prompt and the indirect produced J/tp respectively, which is much more 
involved than that of ip', since we need to consider contributions from the higher charmonium 
states accordingly. And then we shall make a comparison with the newly obtained prompt 
J/ip data by TEVATRON CDF collaboration [31 . 

Assuming the produced J/ ip is measured with the momentum p, we deflne the rest frame 
of J/ip hj a Lorentz boost from its moving frame. The produced J/ip is polarized with the 
polarization vector e* in the rest frame. Within the framework of NRQCD factorization l|, 
after decomposing Dirac- and color- indices, the contribution from the channel through the 
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color-octet (cc)8p5'i] to the differential cross section can be generally written as: 



da[{cc)8[^Si]] = Hij ■ Tij{€, e*, p), 



(1) 



where p = p/|p|, is the 3x3 spin density matrix for producing (cc)8p5'i] and Tij is 
the spin density matrix for the transition of (cc)8p>S'i] into a polarized J/ip, which can be 
decomposed as j?]: 



where aQ i and Co,i,2 are un-determined, non-perturbative but universal parameters. By 
taking the heavy quark spin symmetry, only Cq = {0\Oi^^{^Si)\0) /Q is non-zero, which 
implies that J/ip will have the same spin as the color-octet (cc) pair. In the present, we shall 
make a detailed discussion on the direct J/ip production through pp J/ipln] + X, where 
n stands for the intermediate (cc) -charmonium states up to corrections, i.e. n = (^5*1)1, 
(^5*1)8, (^5'o)8 or (^Pj)s. Further more, the hard subprocess of the hadronic process is 
ab — > J/i'[n] + X, where ab = gg, gq, gq and qq respectively. Based on these processes, 
the relative importance of the undetermined parameters ao,i and Co,i,2 shall be discussed. 
The prompt J/ polarization shall also be discussed, whose signal includes J/ ip meson that 
comes from decays of the higher charmonium states Xa, Xc2 and ip'. Further more, we 
shall take the same spin-flip effect to study the indirect J/ip production from 6-decays, i.e. 
b-^J/ilj + X. 

The paper is organized as follows. In Sec. II, we present the calculation technology for 
the J/ip production. Numerical results and discussions are presented in Sec. Ill, where the 
results for the direct, the prompt and the indirect J /ip polarization shall be presented. The 
flnal section is reserved for a summary. 

II. CALCULATION TECHNOLOGY 

We adopt the same calculation technology as pointed out in Ref . \L\ to calculate the direct 
J / ip production, and for self-consistency, we present the calculation procedure in more detail. 

It is more convenient to transform the spin density matrix T^j into a covariant form. For 
such purpose, we introduce a Lorentz boost matrix Lf , whose components can be written 



Tjj(e, e*, p) = 5ij (e ■ e*ao + e ■ p ■ e* ■ pai) + (ejC* + eje*)co 

+ [{tiPj + tjPi) e* ■ p + (e ^ €*)] ci + piPjE ■ e*C2, 



(2) 
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as [8|J 

4° = f^, (3) 



'0_Pi rj_. I / P^PJ \ Ep 



With the help of such boost matrix, one can transform a purely space-like four-vector, such 
as e = (0, e), from the rest frame of J/ip where the components of p are (M, 0) to the frame 
in which its components are = {Ep,p), Ep = y/M'^ + p^. By applying to both 

sides of Eq.(l2]), and noting the fact that 

g.^Lf^ = -S,„ LfL'f = -g,, + ^, (4) 

we obtain 

Pi^Pp 



Co 



M2 

+ [(cmP^- + (^uP^i) e* ■ p + (e ^ e*)] (-Ci) + {p^^pyt ■ e*) (-C2) H , (5) 



where in the J /ip rest frame, = Lfci = (0, e) and = Lfpi = (0, p). Further more, for 
a particular polarization state A = (0,±1), T^^{e,e* ,p) changes to T^^{e{\),e*{X),p). 

A convenient measure of J/ip polarization is the variable a = (1 — 3.^)/(l + C,), where 
^ = aL/{<yL + o't) , o"L and ax stand for the longitudinal and the transverse components 
of the direct hadronic cross section respectively. As for the indirect J/ip production from 
b decays, ^ = Tl/(Tl + Tt) with and Tt stand for the longitudinal and the transverse 
components of the indirect decay width respectively. 

With the help of Eq.dS]), we can write the total differential cross section for the direct 
J/ip production through the process pp — * J/ip^ln] + X as: 

daxipp J/ip^[n]X) = dXadxbfa/pixa)fb/pixb)dcr^^[n,ab]Tf,^{e{X),e*{X),p), (6) 

ab 

where the differential cross sections for the hard subprocesses, a{ki) +b{k2) — > J / ip[n]{p) + X , 
can be written in the following factorization form, 

da^y[n^ ab] ^ Aab[n]g^u + Bab[n]ki^kiy + Cab[n]k2^^k2v + Dab[n]^^ (7) 

n stands for the (cc)-charmonium state up to corrections, i.e. n = ('^5'i)i, ('^5'i)8, (^5*0)8 
and {^Pj)8 respectively, ki and k2 are the momenta of the initial partons and ab = q q, gq, 
gq, qq. The coefficients ^^^[n], Bab[n\, Cab[n\ and Dab[n] can be read from Refs. jol. [lo| ^. To 



^ We have calculated all these channels and found a good agreement with those in Refs-jol. [lot. 



calculate the longitudinal cross section, we adopt the covariant form of the J / longitudinal 
polarization vector, eL(p)^ = --j=^^=^^^^ ^^m)' where Q = Pp + pp is the sum 



of the initial hadron momenta. In Refs 



111 . Il2| . the fragmentation effect has also been 
resummed to the leading logarithms [a^ lnp^/(2mc)^]" accuracy for the production of the 
color-octet ^5*1 state with the help of the Altarelli-Parisi evolution equation. By taking the 



fragmentation effect into account, the value of a shall be further suppressed [llj. In this 
paper, we shall concentrate our attention on the spin-flip interactions and will not take this 
effect into consideration. 

Further more, the hadronic differential cross section can be simplified so as to obtain J/ip 
Pt distribution: 

where (j[n, a6] = a^y[n,ab]T^y{e,e* ,p), and for a particular polarization state of J/ip, 
do-x, T/,r iAx-\ /■ J J x / \ X / /\ 



dpi 



■{pp-^ J/tlj[n]^X) = Y. I dydXafa/p{Xa) fb/p{x'b) yz ^ Mt y 

ah '' XaXfj[Xa 



where 



:^^^^^^^T,.(6(A),6*(A),p), (9) 



^ Vs XaVs-MTcy ' 



with M|, = yp^ + M^, S = {pp + PpY and y stands for the rapidity of J/ip. 

Secondly, we present the formulae for the prompt J/ip polarization. Theoretical pre- 
dictions of the polarization of prompt J /if) are complicated by the fact that the prompt 
signal includes J / ip mesons that also come from decays of the higher charmonium states Xcj 
(J = 0,1,2), and ip' . The unpolarized differential cross-section daf^""^"^^ "^^^ /dpt is simply 
obtained by adding the unpolarized cross-sections of the various direct-production processes 
multiplied with the appropriate branching fractions. While the longitudinal differential 
cross-section da^"^"^^ '^^^ /dpt is much more involved, which equals 

^ nromptJ/^ ^direct J/^ ^y^j ^^^'^Xcj 

— = — \ \ — H , (11) 

dpt dpt dpt dpt dpt 

where following the discussion in Refs. [l^, Q , we have 
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(13) 
(14) 



and 



+ 
+ 



17 



direct ip' 
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direct ip' 
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-0": 



direct i/)' 
T 



B{^' ^ Xc2 + 7)S(Xc2 ^ W + 7)- (15) 



(5"2,(n) and (3"t(^) stand for the longitudinal and transverse cross sections without the matrix 
elements accordingly, and the summation in the first equation is over n = ('^S'i)i, (^6*1)8, 
(^So)s and (^Pj)8- The direct ijj' production q-'^"^''^^ = Y.n^L{'n){0'^' in)) can be obtained 
from that of J/ip by changing the J/ip matrix elements to ip' matrix elements. The cross- 
section for n = {cc)i{^Pj) is much more involved and we put some necessary formulae in the 
APPENDIX A. 

Finally, by taking the spin-flip effect into consideration, we recalculate the J/ ip production 
from the 6-decay process, h — > J /il)[n\ + X. Following the same procedure of Ref.[l^, the 



total unpolarized cross-section can be written as 



where 



a = (^3(C+ + C^f 

+(2C+ - C-f- 
h = 3(C+ + C-f 



1447r 



cb 



nirmu 1 



4m^ 



X 



a 1 



5m' 



nit 



+ b 



(16) 



[2co + 3ao + ai] + 



2m^ 



Po)) 



2m1 



(17) 
(18) 







^"^3.4 






< 






1.4 ±0.1 

GeV3 


3.9 ±0.7 

lO-^GeV^ 


6.6 ±0.7 

10-2GeV3 


6.7 ±0.7 

lO-iQeV^ 


3.7 ±0.9 

lO-^GeV^ 


0.78 ± 0.36 

10-2GeV3 


9.1 ± 1.3 

10-2GeV5 


1.9 ±0.2 

lO-^GeV^ 



TABLE I: Adopted NRQCD matrix elements from Ref. 
The Wilson coefficients C+(mb) = 0.868 and CJmb) = 1.329 



121. 



151] • In Eq.( fT6il . upon sum- 



ming over the hght quarks s and d, we have used iKbKlP + iKfcKdP ~ iKtP- We have 
taken the initial 6-quark to be unpolarized. Here {O^/^) = {0\O^/^\0) are NRQCD J/ijj 
production matrix elements. If neglecting the spin-flip effects, we return to the same results 
derived by Ref.|l^. As for the longitudinal cross-section from b decays, we obtain 



Ti^ib^ J/ij + X) 



imt 



4m 



2\2 



mlrric 



cb 



mi ± 8m; 



8647r 

± 9(m^ ± 8m^)(C+ + C-f \\ao + ai] ± 



:[2co±4ci±C2] 



mt 



In these equations, we have adopted the relation {0"^'^ i^Pj)) = (2 J ± 1){Oq^^{'^P(^ 



'■ (19) 



III. NUMERICAL RESULTS AND DISCUSSIONS 

As for numerical calculation, we adopt the matrix elements derived by Ref.[12|, which are 
shown in TABJU To be consistent, the parton distribution function is chosen to be CTEQ5L 
[igI and the value of is evaluated from the one-loop formula using the corresponding value 
in CTEQ5L for Aqcd- f^ic = l-5GeV and Mj/^ = 2mc. Both the factorization scale and the 



renormalization scale are taken to be the transverse mass of J/ip, i.e. fif = fir = y M'^ +Pt- 
The collision center of mass (CM.) energy is 1.96 TeV. The branching ratio of J/'ip 
is p = (5.93 ± 0.06) X IQ-^ [17]. As for the charmonium, we take v"^ = 0.30. 

There are two types of power counting rules for the undetermined parameters ao,i and 
Co,l,2 , i.e. the flrst type is 



and the second type is 



oo 2 «i '^1 C2 3 , . 

— ~t> — = — = — ~t> f 20) 

Co Co Co Co 



^^^^^^^^ 0(1). (21) 
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FIG. 1: -R[g5, (^S'i)8] -distributions defined by Eq. (f23]l . where the left diagram is derived by sum- 
ming over all the polarizations and the right diagram is only for the longitudinal polarization of 
J/ip. For the right diagram, the cures for uq and ai are coincide with each other. 

Since the value of v is generally not small for the case of charmonium, those higher v- 
suppressed terms can have a significant impact on theoretical predictions even for the first 
type of power counting rule. 



A. relative importance among the different terms in T^y{e,€* ,p) 

This subsection is served to show the relative importance among the different terms in 
Tfj_i^{e,e* ,p). In order to show clearly the contributions from each parts of T^i/(e, e*,p), we 
rewrite ^^^^ as {n = (^S^s), 

da\ab,n\ da'^°\ab,n\ da'^^lab^n] da'^'^\ab,n\ da°'°\ab,n] da°'^\ab,n] 

, = ^-^Co + j-^ci + j-^C2 + j-^ao + 22 

dpt dpt dpt dpt dpt dpt 

where '^'^ "Jpt'"^^ stands for the pj-distribution for ci = C2 = ao = ai = with an overall 
factor Co being contracted out, and etc.. To show the relative importance among different 
distributions, we define the ratio, 

do-^i [ab,n] da^^ [ab,n\ do-°0 [ab,n] dcr"! [ab,n\ 

dpt dpt dpt dpt 

(23) 

where ab = gg, gq, gq and qq. 
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The i?-distributions for the dominant gluon-gluon fusion mechanism {ab = gg) are shown 
in Fig.([T]), where the left diagram is derived by summing over all the polarizations and the 
right diagram is only for the longitudinal polarization of J/ip- The left diagram of Fig.([T]) 
shows that when summing over all the polarization vectors of J/ip, the weights of and 
Oi shall always at the same order of that of Cq, and more explicitly R'^°[gg, (^5'i)8] = 3/2 
and W^^lgg, (^5'i)8] = 1/2; while those of ci and C2 drop down quickly with the increment 
of pt {0{l/p^)). The right diagram of Fig.([T]) shows that for only the longitudinal part, the 
weights of Oq and ai increase quickly in comparison with that of Cq with the increment of pt 
0{pf); while the weights of Ci and C2 have the same order of that of cq, or more explicitly 
R^'lgg, {'S,)s] = 2 and R'^^lgg, ('S,)s] = 1/2. 



San cii 2ci C2 

IH + — ^ - + — 

2co 2co Sep 2co 



(24) 



B. a simple discussion on the color-octet ^Si matrix element under the spin-flip 
interaction 

Summing over the polarizations on both sides of Eq.(l2]), we obtain a new matrix element 
for the color-octet ^Si state, 

{o\oi/'^Cs,)\oy = 6co 

In principle, the value of the new matrix element {0\Og^'^{^Si)\0)' defined in Eq. (1241) is 
different from that of the usual (0 108^^^(3^1)10), since it involves an extra gauge link in the 
definition and it also takes the spin-flip interaction into account. By taking into the spin-flip 
interaction and the extra gauge links, the matrix element {0\Og^^ {^Si)\0)' contains more 
non-perturbative parameters to be determined, i.e. oq, ai, Cq, Ci and C2. In the following, 
we shall discuss the spin- flip effects based on the two counting rules ( |20l) and (I2T1) . and to 
predict the polarization, we make the ansatz that ai = ci = C2 and introduce two parameters 
/o = flo/co and /i = ai/co, we obtain 

(Oi^'^eS,))' = Qc, (1 + ^/o + . (25) 

fi = "^/o for the power counting rule (12(1 . fi = /o for the power counting rule (1^ respec- 
tively. By setting /o = /i = 0, we return to the result without taking into account the 
spin-flip interaction, which leads to Cq = (0|O8^'^('^S'i)|0)/6. 

In Fig. ([2]), we show the pt-distributions with different value of /o and /i = vfo for the 
direct production of J/ip through (00)8^5*1] at TEVATRON run II {y/s = 1.96TeV), where 




FIG. 2: p(-distributions for the direct production of J /if) through (cc)8[^5'i] at CM. energy ^/s = 
1.96 TeV. Left diagram is derived by summing over all the polarizations and the right diagram is 
only for the longitudinal polarization of J/ip. (3 stands for the branching ratio J/ip ^ f^^l-i- - In 
the calculation, the contributions from all the considered subprocesses have been summed up, and 
rapidity cut \y^/'^ < 0.6| is adopted. 

the left diagram is derived by summing over all the polarizations and the right diagram is 
only for the longitudinal polarization of J/ip. In the calculation, the contributions from all 
the considered subprocesses (with ab = gg, gq, gq, qq) have been summed up, and the rapidity 
cut \y^/''i' < 0.6| is adopted. Fig. ([2]) shows that by considering the spin-symmetry-breaking 
factors ao,i and Ci^2; the pt-distributions derived by summing over all the polarizations of 
J/ip shall not be affected too much in comparison with the case of not considering the spin 
flip interactions (with /o = /i = 0). While for the longitudinal pt-distributions, the spin- 
symmetry-breaking factors ao,i and ci^2 might be important in the large pt regions, i.e. they 
can raise the longitudinal contributions to a certain degree. So such spin symmetry breaking 
corrections can not be ignored safely as has been done in Ref . fiol . since as shown in the 
left diagram of Fig. ([2]), they can change the fraction of the longitudinal part dramatically. 
To show this point more clearly, we present the differential cross-section formulae for J/ij) 
production channel through g + g ^ J /ip{(^Si)^) + (? in APPENDIX B. It can be found that 
Co comes into contributions at 0{l/p1) for the longitudinal differential cross section. This 
is the reason why by taking the spin-symmetry, = cti = Ci = C2 = 0, the longitudinal 
contributions should be neglected at large pt regions, i.e. the J/ip is transverse polarized 
at large pt regions. And it is clear that longitudinal differential cross section will not be 
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suppressed by l/p^, if one takes spin-flip interaction into account, i.e., if those coefficients 
beside cq, especially oq and ai, are not zero. 

In the above, we have shown that the spin-symmetry-breaking factors Oq, cli, C\ and 
C2 might be important in the large pt regions for the production through (cc)8[^>S'i], i.e. 
they can raise the longitudinal contributions to a certain degree. In the large regions, 
it is the (cc)8[^5'i] production channel that yields the transverse polarization, while the 
(cc)8[^5'o] and {cc)^^P.j\ channels both yield unpolarized quarkonia in this limit ^. So, the 
large transverse polarization might be considerably changed by the spin-symmetry-breaking 
factors ao, fli, c\ and C2, and then the large discrepancy between the theoretical prediction 
and the experimental data might be compensated. 



as: 



C. Numerical results of a for direct J j^j production 

The polarization is predicted with the parameter a as a function of j9_l, which is defined 

^_ ( datot _ / f d^tot ^ doA _ 1 - 3^ 

\ dp± dp± J ' \ dp± dp± J 1 + C ' 

where ^ = j^f{j^- If « = 1? the produced is transversely polarized. If a = —1 the 
produced is longitudinally polarized. 

After summing over all the above mentioned channels, we represent the longitudinal 
polarization fraction ^ for the direct J/ ip production as 



Mr. 



^(Oi'*i^S.)) + (»S,)> + + (1 - X) 

^(/o /l) = ^ ' 



(27) 

dpt 

for the case of J/ip through color-singlet (3Sj)-charmonium state, and so on. M:JJ'^ 



where stands for the differential cross section without the matrix element {Oi^^ iP^Si)) 



{Oi'^{^Sa)) + r{Oi'^{^PQ))/ml with r = 3.4 jjj, and the parameter x = {Oi''\^So)) / M'^'^ , 
whose center value is 1/2 and can be varied between to 1. 



^ We have also calculated the spin- flip effects for (cc)8[^S'o], which is quite small in comparison to these of 
(cc)8p5'i], so we shall not take these spin-flip effects into consideration. 
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Pj(GeV) 



Pj(GeV) 



FIG. 3: a for the direct J/'0 production. Left diagram is for f\ = ti/g and the right diagram is for 
/i = /q. The upper shaded band is for /o = 0, the higher middle yellow band is for /o = f^, the 
lower middle green band is for /o = and the lowest band is for /o = 1, where the uncertainties 
of the matrix elements are also included. The upper edge of each band is for x = \ and the lower 
edge for x = 0. The rapidity cut \y'^l'^\ < 0.6 is adopted. The experimental data is from Ref.^]. 

In Fig. ([3]), we present the value of a as a function of pt for the direct J/ip hadronic 
production, where typical values for /o and fi are adopted and the uncertainties from the 
matrix elements are also included. From the figure the produced J/tp will be dominantly 
with transverse polarization at large p±, if one does not take the spin- flip interaction into 
account, i.e., /o = /i = 0. Increasing /q and fi from to 1, a will be decreased accordingly. 
And if one takes /o and /i at the order of 1, a shall be close to 0, which implies that the 
J/ip is unpolarized. 



D. Numerical results of a for the prompt J/ip production 

In calculating a for the prompt J/ip production, we need to know the non-perturbative 
matrix elements {O^"'' (^Pj)) and {O^"' {^Si}) . And for their values we adopt the following 
relations: 

{Of'i^Pj)) = {2J+l){OrCPo)) (28) 

(or^es,)) = {2j+i){oresi)) (29) 

(^^(^Po)) = (9.1 ± 1.3) X 10-2 GeV^ and (Of "(^^o)) = (1.9 ± 0.2) x 10-^ GeV\ As for 
the relevant branching fractions listed in Eqs. (ll3|14|[T5|l . we adopt the values from Ref . jl?! . 
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FIG. 4: a. for the prompt J/i/^ production. Left diagram is for f\ = v/g and the right diagram is 
for /i = /q. The upper shaded band is for /o = 0, the higher middle yellow band is for /o = f^, the 
lower middle green band is for /o = and the lowest band is for /q = 1, where the uncertainties 
of the matrix elements are also included. The upper edge of each band is for x = \ and the lower 
edge for x = 0. The rapidity cut \y'^l'^\ < 0.6 is adopted. The experimental data on the prompt 
J ftp is from Ref.^]. Note the upper shaded band is close to Braaten's results [l^ . 



i.e. Bixco ^ J/^ + 7) = (1-30 ± 0.11)%, Bixa ^ J/^ + 7) = (35.6 ± 1.9)%, i3(xc2 ^ 
J/^ + 7) = (20.2±1.0)%, 13{^' J/^+X) = (56.1±0.9)%, B{^' ^ Xco+l) = (9.2±0.4)%, 
Xci+7) = (8.7±0.4)%, ^ Xc2+7) = (2.6±0.4)%. As for the direct production 
^direct J/-0 ^j^^ ^direct 1/)'^ only the spln-fllp effects in the channel of (cc)8p5'i] are sizable. So 
we shall only consider the spin- flip effects in the direct production cr^iirect j/V' ^^^^ ^dnectip'^ 

In Fig.(j4]), we present the value of a as a function of pt for the prompt J/ip production, 
where typical values for /o and /i are adopted and the uncertainties from the matrix elements 
are also included. The experimental data on the prompt J/ is from the TEVATRON CDF 
collaboration [3]. It can be found that the results for /o = /i = 0, i.e. without taking the 
spin-flipping effects into account, is consistent with Braaten's results 12|]. It is clear that 
under the proper spin-flip interaction, a can be more closer to the experimental data than 
that without these interactions. At large pt, a is reduced by ~ 50% for /g = f ^ and by 
~ 80% for /o = 1. 



13 




FIG. 5: The predicted ob for b J/ip + X as a function of the unknown parameter x. In the left 
diagram fi is taken as vfo, while the right one is for /i = /o ■ 

E. Numerical results of as for indirect J/ip production through b J /if) + X 

By varying the matrix elements within the region of TAB HI we calculate the numerical 
results of versus x for the indirect J/ip production through B J/ip + X, where 
X = {Ol^"^ {}Sq)) / M'l^^ . It is found that the predicted J/i}) polarization parameter 
depends weaker on /q and /i than the case of direct and prompt J/ip production. More 
explicitly under the case of /i = f/o, the range of is shifted from [—0.100,0.193] to 
[—0.123, 0.175] by varying /o from to 1, while under the case of /i = /o, the range of as is 
shifted from [—0.100, 0.193] to [—0.143, 0.163] by varying /q from to 1. On the other hand, 
it dependents heavily on the matrix elements {O'l^^ (}Sq)) and M^/^ . To be consistent with 
as derived in literature, e.g. as = —0.13 ± 0.01 for J/ip events with pxiJ/"^) > ^GeV by 
CDF group [20I. a larger x that approaches 1 should be taken, as is implicitly adopted by 



Ref. 



14| . Since the spin-flip effect is weaker in this indirect J/ip production, a more precise 



measurement of it can predict more precise matrix elements. As a cross check, we have found 
that without considering the spin-flip effect and by varying the involved matrix elements 
within the same uncertainty region derived by Ref. [l^ . we can obtain the same allowable 
region for a^, i.e. as € [—0.33,0.05]. 
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IV. SUMMARY 



It is noted that the T{nS) production is somewhat different from the case of J/ip and ip'. 



Ref. 



231] shows that the NLO correction plus the LO resuhs for the color-singlet of T{nS) can 



explain both the total unpolarized and the polarized T{nS) production cross sections. Then 
there is no need to consider the contributions from the color-octet transitions for T{nS) 
production, or in another words, the color-octet components give negligible contributions 
to the T{nS) production. For the J/ip or ip' production, the NLO correction plus the LO 
results for the color-singlet production can also lead to longitudinal polarized J/tp or ip' 
[4]. However, it is well-known that with the color-singlet contribution only, one can not 
explain the unpolarized J/tp or ip' production cross section. And by taking the color-octet 
cc components into consideration, one can well explain the total unpolarized cross section 
of J/ip or ip' 2J], which is regarded as a great triumph of NRQCD. Within the NRQCD 
framework, if keeping the spin symmetry for the charm quark, Refs.0, ^ show that the 
J/ip polarization puzzle can not be solved even by including the NLO corrections. Hence to 
well explain both the unpolarized and longitudinal cross sections of J/ tp or tp' is much more 
involved than the case of T{nS). 

We have shown that the spin-flip interaction can have a significant impact on the transi- 
tion of a color-octet ^Si cc pair into J/ip. Such impact can be parameterized by introducing 
new parameters in the transition matrix T, i.e. a^ i and Co,i,2- If the heavy quark spin sym- 
metry holds, the matrix has only one parameter cq. The newly introduced parameters ao,i 
and Ci^2 are power suppressed in comparison to that of Cq in principle. However the charm 
quark is not heavy enough, or v is not small enough, these new parameters are not small 
in comparison with cq, or even they can be at the same size of cq. Numerically, it is found 
that these parameters can significantly reduce the polarization parameter a of J/ip. More 
explicitly, we have calculated the direct J/tp polarization, the prompt J/tp polarization and 
the indirect J/ tp polarization from the b decays. 

Without the spin-flip effect, we return to the same results with those in literature under 
the same parameters. While by taking the spin-flip interaction into consideration, the pre- 
dicted a is more close to those measured at TEVATRON. At large pt, a for the prompt J/ip 
is reduced by ~ 50% for /o = and by ~ 80% for /o = 1. Then such spin-flip interaction 
as have been argued by several authors may provide a suitable way to solve the J/ip po- 
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larization puzzle at TEVATRON. Since the NLO correction shall provide a large K factor 
of the total cross section (ratio of NLO to LOV, e^. i^' ~ 2 for color-singlet (cc)i[^5'i] and 
~ 1 for color-octet (cc)8[^5'o] and (cc)8p5'i] [4, l5|, and because the NLO can increase the 
transverse distributions of produced J/ip more than the total distributions, then the value 
of a can be further lowed by including NLO results into our present calculation. Further 
more, it has been argued that the production of J / ip associated with a cc quark pair might 



also help to dilute the J /if) polarization 22j. Such analysis is out of the range of the present 
paper, which is much more involved since it involves a NLO calculation of these processes 
with spin-flip effects being under consideration and a newly systematical determination of 
the color-octet matrix elements. More over, we have found that the predicted indirect J /if) 
polarization parameter ub depends weaker on the spin-flip effects than the case of direct and 
prompt J / ip production. And then a more precise measure of the indirect J / ip polarization 
from the h decays can be adopted to predict more precise color-octet matrix elements, which 
can inversely improve our estimations on the direct and prompt J/ip production. 
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V. EXPAND '^^'?^y9a] rpjj^ LARGE pt LIMIT 

A. Formulae for the production of J/ip through the channel of (cc)ipP2] 
For the case oi n = (cc)ipP2], the differential cross-section takes the form: 

daxipp^ J/i'^in)X) = dxadxbfa/p{xa)fh/p{.Xh)da^,upa[n,ah]p^^^^", (30) 

ah 

where we have 

'^^^''^d!^' '^^^ ^ ^afeh]fi'/xpfi'i/a + Bab[n]g^pki,yki^ + C ab[n]g ^,pk2yk2a + Dab[n]g^pki,yk2a 
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where 



and 



and 



+Eab[n]ki^kipki^ki„ + Fah[n]k2pk2pk2yk2a + G ab[n\kif,kipkiyk2a 

+Hab[n]k2pk2pk2ukla + Iab[n]kipkipk2uk2a + Jab[n]kipk2pki^k2a, (31) 



2 1 1 

A=~2 2 3 



= e^-*(0)e^'^(0) = i(2pr " Pr)(2pr - pD (33) 

pr" = E e'^^*(A)e^-(A) = ^(p{;vr + prpr + Po>r + poVD (34) 

|A|=1 ^ 

pr"^ = E e'^^*(A)e'''^(A) = lip^pT + pTpT - P^'^pT) (35) 

|A|=2 ^ 



= ^ 6'^*(A)6'^(A) = -(^^'^ + ^ (36) 
p- - 'z^l , e^(0) ^ ^ (37) 

/(p • QY - M2g2 



pr = Ee"*(A)e^(A) = p'''^-pr- (38) 

fci and k2 are the momenta of the initial state partons, ah = gg, gq, gq and qq. P and Q 
are the momenta of the bound state and the total four-momentum of the colliding hadrons 



respectively. All coefficients, ^^^[n], Bab[n], Cab[n] and etc. can be read from Refs.|9l. Illl. |21|. 



B. Formulae for the pt-expansion of the dominant gluon-gluon fusion subprocess 

As for the dominant subprocess: g{pi)+g{p2) 'J/i'{{^Si)s){p3)+g{pA), in the laboratory- 
Frame, we have 

Pi = —ixa,0,0,Xa), P2 = — (a;fc, 0, 0, -a;^) , p^ = (M^ cosh(?/),p^,p|, sinh(y)), 

(39) 

where pi = {p'i , pf , pf , Pi) , y is the rapidity of J/ip, = +Pt, and we have 

S = (p^+p^f = XaXbS, 

i = (Pi - Ps)' = - VSMTXa[cosh{y) - sinh(y)], 

u = (p^ - p^y = - v^MTX6[cosh(?/) + sinh(2/)], (40) 
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with S the square of CM. energy for the hadronic coUider. Using the above formulae, we 
obtain the pt-expansion for the dominant gluon-gluon fusion subprocess in the large pt limit: 



dt 



1 

Pt 
1 

pi" 



3eyh7iW{M\xa - e^yxhf - SxaXb{xa + e^yx.Y) 



2M^Sixa^Xb^{xa + e^yxby 



36(1 + e^y?M^S^xJxbHxa + e^yxbf 



27(/2 - fl)S'XaV - 



'27e'^Hfi - f2)S^Xa^Xb'' + 2e'''yxbHfiXail00M^ - SlM^SxaX^ + 

27S^Xai2Xa - Xb)Xb^) + 27 f^Sxh" {-{M^ X^) + SXa\2Xa + Xb))) - 
2e^yXa^{flXb{-100M^ + SlM^SXaXb + 27S^Xa^{Xa - 2Xb)Xb) + 

27 f2Sxa\M^Xa " Sxb\xa + 2xb))) + 6^^'^' (/i (-8M^ (Slx^ - 50a;,,) 

-?,2AM^SXaXb{Xa + Xb) + 27 XaXb^ {IQXa^ + Sx^Xb - Xb^)) + 

27 f2Sxb''{-4:M^Xb + SxaiGxa" + SxaXb + Xb''))) + e^yx^^Xb ■ 

{fi{8M^{50Xa - SlXb) - 324:M^SXaXb{Xa + Xb) - 27S^Xa^Xb{Xa^ - 
8XaXb - 10 V)) + 27f2SxJ{SXb{Xa'' + SXaXb + 6Xb^) - 4:M^Xa)) 
+2e^yXaXb(-27f2SXaXb(M^{xJ + Xb'^) - 2SXaXb(xJ + SXaXb + Xb')) 
+fl{-81M^SXaXb{Xa' + AXaXb + Xb') + BAS^xJxb'ixJ + 5XaXb + Xb') 



+4M\25xa' - 62xaXb + 25xb^))) 



+ 



.pI 



(41) 



where /i = 3ao + ai + 2co and /2 = 4ci + 3c2. It can be found that /2 comes into contribution 
at least at 0{l/pf) in comparison to /i. While for longitudinal distribution of J/ip, we only 
need to make the change: fi gi {i — 1, 2) with gi — ao + ai, g2 — 2cq + 4ci + C2, i.e. 

daL[{'S,)s,gg] da[^Sl'\gg] 



ifi 91, f 2 g2)- 



(42) 



di di 

One may also find that for the longitudinal part, Cq comes into contributions at 0{l/p^) in 
comparison to the total summed results, this is the reason why by taking the spin-symmetry, 
i.e. Qq — Qi — Ci — C2 — 0, the longitudinal contributions should be neglected at large pt 
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regions, i.e. the J/ if) is transverse polarized at large pt regions. 
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